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In this Letter we study the effect of Pauli blocking on Efimov states in a quantum Fermi gas and 
illustrate that the universal Efimov potential is altered at large distances. We obtain the universal 
spectrum flow of Efimov trimers when the Fermi density is varied and further consider the effect of 
scattering of trimers by the Fermi sea. We argue that the universal flow is robust against fluctuating 
particle-hole pairs that result in an infrared catastrophe in impurity problems. 



Recently, universal three-body Efimov structures that 
were first proposed in the 1970s have been successfully 
probed and determined in a series of remarkable cold- 
atom experiments on inelastic loss spectra (H-[l3j. The 
current data that cover an impressive range of scatter- 
ing lengths quite conclusively demonstrate the univer- 
sality of the Efimov structures [TH, 0]. The spectacu- 
lar progress made on this subject has further illustrated 
that the physics of cold atoms at sub-microkelvin tem- 
peratures can shed light on fundamental issues in other 
systems at very different energy scales, in this case the 
few-body structures of nuclei. 

Although Efimov physics 0, 0] and the theory of loss 
rates (3, |4j have so far been quite successful in explaining 
many aspects of the data, in experiments on cold gases 
there are always many-body backgrounds. Logically, it 
is important to understand how few-body states respond 
to the presence of a quantum many-body background be- 
low degeneracy temperatures. For instance, can Efimov 
trimers survive a many-body background and how are 
they affected when scattered by the background? An 
equally important and challenging question is what kind 
of many-body correlations can be induced by the uni- 
versal few-body states in a quantum gas or mixture. In 
this Letter, we attempt to answer one of these questions, 
and particularly to examine the effect of Pauli blocking 
of scattering or open-channel Fermi atoms on the three- 
body Efimov states. To answer this question, we solve 
the three-body spectrum in the presence of a Fermi sea 
and focus on the effect of Pauli blocking on the Efimov 
physics. These spectra are known to be one of the cor- 
nerstones for the theory of three-body recombination and 
the loss spectrum of metastable quantum gases. Our re- 
sults can also be applied to study the dynamics during 
the initial stage of recombination and estimate the life- 
time of quantum gases. 

In the following, we consider that two identical heavy 
Bose atoms with mass M have an interspecies resonance 
with a light Fermi atom with mass m(<C M) in the 
presence of a Fermi sea of light atoms. In the Born- 
Oppenheimer approximation, we can analyze the fast mo- 
tion of the light atom assuming that the two slow heavy 
atoms are a distance R apart. In the absence of a Fermi 
sea, simultaneous near-resonance scattering of the light 



atom by the two heavy atoms induces a bound state, and 
the binding energy at resonance is equal to h 2 fl 2 /2fiR 2 , 
where fl = 0.567 is the root of O = e~ n (see below) and 
fi = mM /(m + M) is the reduced mass of a light-heavy 
subsystem. The bound state therefore glues together the 
two heavy atoms and yields an attractive long range po- 
tential — h 2 tt 2 /2/iR 2 . This universal l/R 2 potential plays 
a paramount role in Efimov physics and results in the 
spectacular universal hierarchy structures in three-body 
spectra P, Q • We then investigate the effective attractive 
potential between two heavy atoms in the presence of a 
Fermi sea with Fermi momentum Tik-p. The interspecies 
interactions can be treated as zero-range ones [3, EH , so 
that for a light atom interacting with two heavy atoms 
at ±R/2, V(r) = V [8(r - R/2) + 6(r + R/2)] where 
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EqfT] relates the contact interaction strength Vb to the 
scattering length ohl via a standard regularization proce- 
dure [16|. Here e£ = h 2 k 2 /2fi and now the reduced mass 
fj, ~ m since M is much heavier than m. The Schrodinger 
equation for </>(k) , the momentum-space wavefunction for 
the light atom, is 
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is the kinetic energy for the light atom, 



and (/>(k) = when k is within a spherical Fermi surface 
of radius fcp due to the Pauli blocking effect. We find 
that the binding energy is given by 
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where we have set i*l = 2raE\ J /h 2 . This equation is 
valid for — oo < i*l < k F , and has a unique solution for 
arbitrary ohl- One can easily verify that the solution 
can be written as ul(«hl, -R, &f) = — fcp5(^FaHL, kpR), 
where g is a dimensionless scaling function. 
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FIG. 1. Effective potentials u^/kp versus k-pR. From top to 
bottom we plot the potential for /cfcihl = — 1, oo, 1. At reso- 
nance, Uh(R) = —Q 2 /R 2 if the Pauli blocking effect is absent 
(dashed line); the dotted line is the estimated potential due 
to the scattering of trimers by the Fermi sea (see discussions 
on page 4). 



In the absence of a Fermi sea or when fcp — > 0, ul 
satisfies the following simple equation: y/—ui~ = + 

e v/ ff " L ■ So for a positive scattering length and R = oo, 



El = — h /2mog L and the light atom forms a two-body 
bound state with one of the heavy atoms. However, when 
R <C I^hlI the bound state is severely affected by the sec- 
ond heavy atom and El = — h 2 fl 2 /2mR 2 , even when ohl 
is negative. One can also show that right at resonance, as 
a result of non-linear interference between waves coming 
off the two heavy atoms, the effective scattering length 
for the light atom is proportional to R, the distance be- 
tween the two heavy atoms. 

Now we turn to the effect of a finite density of fermions 
and focus on the resonant case where 1/ohl = 0; in this 
limit g is a function of fcpi? only. At short distances when 
kvR <C 1, the effective attractive potential ul = — ri 2 / i? 2 
and is determined by the motion of a single light atom. 
However, this universal behavior is completely changed 
when kpR becomes order of unity or larger and the col- 
lective effect becomes important. First at kpR — 0.799, 
the bound state energy El increases to zero. After this 
point, El continues to rise until it eventually reaches a 
maximum and then decreases, settling into a pattern of 
oscillations around a saturation value; its asymptotic be- 
havior is 
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where s» 0.695fcp. See FigJT] where details are shown 
numerically. Note that h 2 u L °/2(j, is identical to the two- 
body binding energy at resonance and therefore repre- 
sents the atom-dimer threshold for trimers. 

To understand the Pauli blocking effect on trimers, we 
adopt a k-space approach instead of employing the usual 



hyperspherical coordinates [l[ . The Schrodinger equation 
for a three-body wavefunction <I>(ki, k 2 , k 3 ) is 
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where ki is the momentum of the light atom, k 2 and 
k 3 are the momenta for the heavy atoms and e^ 1 = 
h 2 kl/2m, = h 2 \4/2M. $(k 1) k 2 ,k 3 ) vanishes when 
ki is within the spherical Fermi surface of radius fcp ■ Fur- 
thermore, $(ki, k 2 , k 3 ) = <£>(ki, k 3 , k 2 ) because of the 
Bose atom exchange statistics. For now we will consider 
the zero-total-momentum subspace and further introduce 
a wavefunction for one heavy atom relative to the heavy- 
light dimer formed by the other two atoms: 

*AD(k 2 ) = / dq $(q,k 2 - q, -k 2 ). 

J\q\>k F 

The resulting function ^AD(k) is shown to satisfy 
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where 7 = M(M + m)/(2M + m) is the reduced mass for 
collisions between a heavy atom and a heavy-light dimer; 
also e 7 — '- 2 
as 



k = h k 2 /27. Integrals J are over a region defined 
w k + q| > fcp to exclude the occupied states. 
When ^ approaches 1, or in the Born-Oppenheimer ap- 
proximation, we can integrate out the fast degrees of the 
light atom and map the three-body problem to a simple 
equation for the dimer and heavy atom: 
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Eq(7]is valid for —00 < u B < fcp- Here u B = 2fiE/h 2 
and (3 — y/j/fi which is much larger than one in our 
case. The kinetic term (the square-root term on the left 
hand side of the equation) scales linearly as a function 
of k for k ^S> /3y/uB, reflecting the composite nature of 
the dimer at short distance; this peculiar structure was 
appreciated in earlier studies on dimer-atom scattering 
14 Il7| . One can also write down the corresponding dif- 
ferential equation for the Fourier transformed wavefunc- 
tion 'F ad (r) = / <ik<FAD(k) exp(ik-r). In the limit where 
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FIG. 2. (color online) A subset of eigenvalue trajectories at 

resonance (ahl = oo) in the Z — sgn(uB)y/\uB\ro versus 
X = fcpro plane; the red trajectories are for trimers dressed in 
fluctuating particle-hole pairs estimated in the leading order. 
The ray Z = \/Y°°X (dashed line) represents the dimer-atom 
threshold. The termination points along this threshold X!^ 
are discussed in Eq[TD] Here ft = \p~ijyi = 10, and ro is the 
hard-core size of Bose atoms. Only Efimov states are shown, 
and we take 5 — n/A [T^ ]. 



f3 is infinite, Eqj7]for V^ad^ = R-) is equivalent to Eq.3 
for two heavy atoms at a fixed distance |R| = R apart. 

To obtain the spectrum flow of Efimov states, we carry 
out a WKB analysis of the wavefunction f AD(r) = 
exp(i/3So(r) + iSi(r)) and the bound state energies 
Mb; the WKB approach is valid as far as j3 (~ 
\/M/2m) is much larger than unity. By comparing 
the resulting equation for So, derived from Eq[7l to 
EqlHl one further establishes that So is a simple func- 
tion of Ui,(r). For the S'-wave channel, "JadM = 
e ' Sl r ir) sin (ji J^dr'y/iiB - u^(r) + S^j for r < r\, where 
r\ is the first semiclassical turning point defined by 
ul{ti) = ub [3 and 6 is a phase shift calculated be- 
low. When r approaches zero, the phase integral has a 
logarithmic divergence for any value of k-p, as seen in 
the hyperspherical approach to the three-atom problem. 
Without losing generality, we proceed by introducing a 
Bethe-Peierls boundary condition at r = ro to take into 
account a hard-core repulsion with range ro that is much 
shorter than the Fermi wavelength 2ir/k-p. The quanti- 
zation condition is: 

u f ri 1 

w(k F r 07 tj) = / dr'y/u a - u L (r') = -(mr - S) (8) 

where n = 1,2,3,... indexes the eigenvalues. Since 
Uh{r)/ kp is only a function of kpr, the dimensionles phase 
w is only a function of kpro and mb/^f- ^ can be ob- 
tained by applying the standard matching formulas for 
the Schrodinger equation, and if there is a single turning 
point we get 6 = 7r/4. fig ]. 

The simplest case is where fcp = 0. Then there are 



infinitely many solutions that accumulate near zero en- 
ergy, and u n /un+i — > as n — > oo, with eS w 254.5 
[110 ■ F° r a non-zero fermion density, we introduce three 
dimensionless parameters: X — k-pro, Y = ub/^f ana - 
Z = sign(Y)y / \Y\X . Our results on the spectrum flow 
are presented in FigfJ] in the Z-X plane, where each 
eigenvalue trajectory shows how the solution to Eq.8 for 
fixed n varies as we increase X. For fixed non-zero X, the 
ratio u n /u n+ i depends on n and approaches the univer- 
sal value of exp(27r/f2/3) only for low lying states with n 
much less than N(X), where N(X) denotes the number 
of Efimov states at a given density X. Since N(X) is a 
rapidly decreasing function of X , the universal value can 
only be attained when X is small. As X increases, each 
eigenvalue increases towards u^, the dimer-atom thresh- 
old, eventually reaching this value and disappearing from 
the spectrum; the eigenvalue trajectories terminate along 
the ray Z = X^u^ /fcp in the Z-X plane. The termi- 
nation point at which u n reaches the dimer-atom 
threshold is given by the solution to: 

w(X™,Y°°) = ~(nn-6°°) (9) 

where Y°° — u^/kp — 0.695 and 5°° is the phase shift 
at Y = Y°°. 

As far as X is much less than unity, the termination 
points have certain universal properties. Indeed, as X 
approaches 0, w{X, Y) approaches an asymptotic form 
w(X, Y) — > — illogX + v(Y) where v(Y) is a regular 
function of Y for — oo < Y < Y°°. This formula yields 

estimates for several quantities of interest. For example, 

»(y°°) 

from Eq.8 we find that N(X) -> SD. \ og e ~^~ + £1 as X 
approaches 0. Furthermore, the termination points are 
X™ ps e n(«( i ' 00 )+ i f ) e ~m for small A^°, which leads to 
a universal relation between different termination points: 

*£W+i = e*. (10) 

for small X^> or large n, and otherwise independent of 
n. We can also calculate X n (Y), the point at which u n 
increases to Y, for — oo < Y < Y°°] we obtain a re- 
sult identical to Eq.10 after replacing X%° / ' X^ +1 with 
X n (Y)/X n+ i(Y) in that expression. These are the new 
universal properties in three-body physics when there is 
a Fermi sea quantum background. Our numerical results 
agree well with the WKB results: for /3 = 10, the differ- 
ence is within a few percent for the densities shown in 

Fig. m 

These universal structures in the spectrum flow should 
be robust against the scattering of trimers by the Fermi 
sea. So far we have treated the Fermi sea as static or 
incompressible, but a trimer can further collide with the 
Fermi sea, exciting particle-hole pairs near the Fermi sur- 
face and polarizing the background. One of the main 
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effects of these fluctuating pairs is to suppress the inter- 
species interactions, as occurs in an interacting Fermi gas 
[20( 1 , leading to a reduction of the Fermi-Bose dimer bind- 
ing energy (estimated in Ref. 2l|). Following those dia- 
grammatic calculations, in the limit of heavy Bose atoms 
one finds that the dimer binding energy measured from 
the Fermi energy is reduced by a factor of m/M when 
/cfOhl is small and negative. This can also be attributed 
to an effect related to Anderson's infrared catastrophe 
22]. As a result, Y°° , which specfies the dimer-atom 
threshold for trimers, moves closer to unity (see FigJIJ. 
Y°° and v(Y°°) are now functions of j3, and these mass- 
ratio-dependent corrections to the effective potential shift 
the flow of each eigenvalue trajectory. However, they do 
not affect the universal relation between different tra- 
jectories, as EqJTU] is independent of Y°° and v(Y°°); 
we obtain the same universal value for X^°/X^_ 1 , pro- 
vided is small. At densities close to and related 
by Eq llOl we anticipate distinct peak-dip signatures in 
the inelastic loss spectrum similar to those observed in 
Ref.liO]. 

In conclusion, we have shown that the spectrum flow of 
Efimov states when the density of fermions varies is uni- 
versal as far as the hard-core size of heavy Bose atoms ro 
is much shorter than the average interatomic distance. 
Our results can be applied to understand three-body 
states in 87 Rb- 6 Li and 23 Na- 6 Li mixtures. It is possi- 
ble to generalize the above idea to other situations. For 
instance, when applying the same approach to the case of 
one light boson in resonance with two heavy Fermi atoms 
in the presence of a Fermi sea, we find that because of the 
anti-screening effect of background Fermi atoms, the me- 
diated potential ul{t) = — fl 2 /r 2 [l — qorsi(qor)], where 
<7o = /3 2 fcp/(27r) 2 and si(x) is a sine integral which os- 
cillates and approaches — — when x — > oo. A 
similar spectrum flow can be obtained and details will 
be presented elsewhere. There are also open questions 
that need to be answered in the future. One is the role 
of hole-like configurations and the effect of Fermi sea 
polarization on the three-body spectrum 23]. Another 
question is what are the dynamic consequences of these 
Efimov states embedded in a quantum gas and how do 
they affect the inelastic loss spectrum and dimer-atom 
or dimer-dimer elastic scattering 17|? The eventual fate 
of such a quantum gas after a large fraction of scatter- 
ing atoms fall into the Efimov channel remains unknown. 
A straightforward way to proceed is perhaps to include 
more atoms in the Efimov channel and extend the analy- 
sis to 4-body, 5-body etc Along that road, quantum 
Monte Carlo simulations have been performed to under- 
stand cluster structures of bosonic atoms 2J, |25j . Sim- 
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